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The observed baryon asymmetry of the Universe can be due to the B — L violating decay of 
heavy right handed (s)neutrinos. The amount of the asymmetry depends crucially on their number 
density. If the (s)neutrinos are generated thermally, in supersymmetric models there is limited 
parameter space leading to enough baryons. For this reason, several alternative mechanisms have 
been proposed. We discuss the nonperturbative production of sneutrino quanta by a direct coupling 
to the inflaton. This production dominates over the corresponding creation of neutrinos, and it can 
easily (i.e. even for a rather small inflaton-sneutrino coupling) lead to a sufficient baryon asymmetry. 
We then study the amplification of MSSM degrees of freedom, via their coupling to the sneutrinos, 
during the rescattering phase which follows the nonperturbative production. This process, which 
mainly influences the (MSSM) _D— flat directions, is very efficient as long as the sneutrino quanta are 
in the relativistic regime. The rapid amplification of the light degrees of freedom may potentially 
lead to a gravitino problem. We estimate the gravitino production by means of a perturbative 
calculation, discussing the regime in which we expect it to be reliable. 
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I. INTRODUCTION AND SUMMARY 

The generation of the Baryon Asymmetry of the Uni- 
verse (BAU) represents one of the puzzles of Cosmol- 
ogy. Three ingredients are required to achieve this 
task: baryon number violation, C and CP violation, and 
departure from thermal equilibrium. The baryon num- 
ber violation can be challenging to implement, because it 
must be consistent with the current lower bound on the 
proton lifetime, Tp > lO-^^ years The Standard Model 
(SM) is a C and CP violating theory, and contains non- 
erturbative B + L violating interactions (sphalerons) 
-which are rapid in the early Universe but unable to 
mediate proton decay. However, it seems difficult to use 
this baryon number violation to create the asymmetry in 
the SM |1 and its more popular extensions An at- 
tractive alternative is to generate a lepton asymmetry |^ 
in some C, CP and lepton number (L)-violating out-of- 
equilibrium interaction, and then allow the sphalerons to 
reprocess part of it into a baryon asymmetry. An appeal- 
ing feature of this scenario is that while neutrino masses 
are experimentally observed (and are L-violating, if 
they are Majorana), there is still no evidence for baryon 
number violation. 

The above idea is naturally implemented in the con- 
text of the see-saw |^ , which is a minimal mechanism for 
generating neutrino masses much smaller than the ones 
of the charged leptons. Three right-handed (r.h.) neu- 
trinos Ni are added to the SM particle content, given 
Yukawa interactions with the lepton and Higgs doublets, 
and large Majorana masses. This gives the three light 
neutrinos very small masses, due to their small mixing 
with the heavy r.h. neutrinos through the Dirac mass. 
Grand Unified Theories (GUT) and their supersymmet- 



ric versions, that constitute natural candidates for the 
Physics beyond the SM, often contain r.h. neutrinos in 
their particle content. In this paper we consider the su- 
persymmetric version of the see-saw mechanism, which is 
theoretically attractive because it addresses the hierarchy 
between the Higgs and r.h. neutrino masses. 

The r.h. (s)neutrinos of the see-saw can generate the 
BAU via leptogenesis, in a three steps process j^, |l^, 
[Tl] , p^ . First, some {CP symmetric) number density 
of (s)neutrinos is created in the early Universe. Then, 
a lepton asymmetry is generated in their CP violating 
out-of-equilibrium decay. Finally, the lepton asymmetry 
is partially reprocessed into a baryon one by the B + L 
violating interactions, provided it is not washed out by 
lepton number violating scatterings. In this paper, we are 
mostly interested in the first step, although in the next 
section we will also briefly review the decay and washout 
processes. 

The most straightforward and cosmological model- 
independent mechanism to generate r.h. (s)neutrinos is 
via scattering in the thermal bath However, as dis- 
cussed in section the parameter space available is re- 
stricted in supergravity-motivated models. Indeed, un- 
less some enhancement of the CP asymmetry character- 
izing the r.h. (s)neutrino decay is present (which occurs 
for example if they are nearly degenerate in mass) the 
generation of a sufficient lepton number poses a rather 
strong lower bound on their mass [|l3|, Q |l^ (see also 
Jl6t). For thermal production, this translates into a lower 
bound on the reheating temperature Tjih of the thermal 
bath. On the other hand, if supergravity is assumed, Thh 
cannot be taken arbitrarily large without leading to an 
overproduction of gravitinos IitI Ii|, |2^, |2|, p|. The 



two requirements are compatible only provided a nearly 
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maximal CP asymmetry (again, banning any possible 
enhancement from mass degeneracy) is present in the 
r.h. (s)neutrino decay. If this is not the case, alterna- 
tive mechanisms of production for the r.h. (s)neutrino 
have to be considered. 

As remarked in leptogenesis can be achieved if 
at least one r.h. sneutrino has a smaller mass than the 
Hubble parameter, i.e. [|7^ M]s < H, during inflation. 
In this case, quantum fluctuations of this sneutrino com- 
ponent are produced during the inflationary expansion, 
and amplified to generate a classical condensate. The de- 
cay of the condensate eventually generates the required 
lepton asymmetry. The above requirement Mjv < H is 
not trivially satisfied in a supergravity context, since su- 
pergravity corrections typically provide a mass precisely 
of order H to any scalar field of the model |Q . In this 
case, however, a suitable choice of the Kahler potential 
can induce a negative mass term m,?^^ ~ —H^, so that 
a large expectation value will be generated for the sneu- 
trino component during infiation . This also leads to 
the formation of a condensate during infiation, and to 
successful leptogenesis as in the previous case. 

Large variances can be produced during infiation if the 
sneutrinos are not too strongly coupled to the infiaton 
field since this would generate a high effective mass 
which could fix {N) — during inflation. However, if 
one of the r.h. (s)neutrinos is coupled to the inflaton, 
there is the obvious possibility that a sufficient amount 
of (s)neutrino quanta is generated when the inflaton de- 
cays. Quite remarkably, for a rather wide range of models 
this decay occurs in a nonperturbative wa y ||25| , (this 
is known in the literature as preheating |p6|). In mod- 
els of chaotic inflation this is due to the coherent 
oscillations of the infiaton field, which can be responsi- 
ble for a parametric amplification of the bosonic fields 
to which the infiaton is coupled ||7^. It is important to 
remark that this resonant amplification docs not require 
very high couplings between the infiaton and the pro- 
duced fields. For a coupling of the form {g^ /2) 0-^ N'^ in 
the scalar potential, resonant amplification of the field N 
already occurs for g^ > 10~® , if the mass of N is neg- 
ligible at the end of inflation, and if a massive infiaton is 
considered. For a massless inflaton {V {(f) = A(/)'*/4), an 
efficient resonance is present also for much smaller values 
of g (we will show this explicitly in section 3), since in 
this case the resonance is not halted by the expansion of 
the Universe [||. 

If the produced particle is very massive {M^ ft "^0)1 
the effectiveness of the resonance becomes a highly model 
dependent issue. A potential of the form V {(/>) + 
Mfj -I- 5^ (j)^ iV^/2 , has been considered in the liter- 
ature mainly to discuss the production of heavy bosons 
needed for GUT baryogenesis |Q. Working in the 
Hartree approximation, it has been found pl[ that a res- 
onance is effective only provided the coupling g^ satisfies 
g^ > 10"'' (MAf/m^)^ . Taking into account all the other 
backreaction effects, a stronger lower bound on g has to 
be expected |3^ , since the latter typically limits the 



growth of the fluctuations amplified by the resonance. 

Very different bounds can be expected for dif- 
ferent potentials. Consider for example V {(f)) + 
{Mj^ + g (/)) In this case, due to the high ini- 

tial amplitude of the inflaton oscillations, the total mass 
of TV can vanish at some discrete points even for a cou- 
pling as small as g'^ ~ 10"^" [MM/m^f' . Whenever 
Mat + g<i> = , parametric amplification of N occurs. 
Thus, the lower bound valid for the previous potential is 
considerably weakened. Although this second choice of 
the potential may seem ad hoc, we note that it is the one 
which arises in supersymmetric models if both the r.h. 
sneutrino mass and interaction with the infiaton are en- 
coded in the superpotential, W (N) D N'^ + g4> ■ 
We regard this as a very natural possibility. 

The idea of a nonperturbative production associated 
to the vanishing of the total mass has been applied to 
leptogenesis in Q. The analysis of Q focused on the 
production of r.h. neutrinos, with a mass term of the 
form (Mat -|- g (j)) NN . From the results of |Q , and from 
the analytical computations of , it can be shown that 
a sufficient lepton asymmetry is generated if the mass of 
the r.h. neutrinos is higher than about 10^"' GeV, and if 
their coupling to the inflaton satisfies g > 0.03 (we will 
derive these bounds in section 2.3). Here we note that 
this high coupling can in principle destabilize through 
quantum effects the required flatness of the inflaton po- 
tential. This, in addition to the strong hierarchy between 
the r.h. neutrino mass and the electroweak scale, moti- 
vates the study of the supersymmetrized version of the 
mechanism proposed in . 

One of our aims is to show explicitly that, in the su- 
persymmetrized version of the above model, the nonper- 
turbative production of the r.h. sneutrinos is much more 
efficient than the one of the neutrinos. Due to supersym- 
metry, the infiaton couples with the same strength both 
to the r.h. neutrinos and to the sneutrinos, so that if the 
former are produced at preheating this will also occur for 
the latter. However, while production of fermions is lim- 
ited by Pauli blocking, the production of scalar particles 
at preheating is characterized by very large occupation 
numbers. This high production has typically a big im- 
pact on the dynamics of the inflaton fleld. The most 
immediate backreaction effect is the generation of an ef- 
fective potential for the zero mode of the inflaton. This 
effective potential, taken into account in the Hartree ap- 
proximation p^ , is typically comparable with or even 
dominant over the tree level potential V (0) . There are 
however two equally important effects which are beyond 
the Hartree approximation. The first is due to the scat- 
terings of the produced quanta against the zero mode of 
the inflaton. This destroys the coherence of the oscilla- 
tions, thus ending the resonant production characterizing 
the early stage of preheating |Q . The second is the am- 
plification of all the other fields to which the produced 
quanta are coupled. This is a very turbulent process, 
dominated by the nonlinear effects caused by the very 
high occupation numbers of the fields involved. As a 
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result, all these mutually interacting fields are left with 
highly excited spectra far from thermal equilibrium JM . 
Both these effects are denoted as rescattering ||3^ . [[76t 

Rescattering strongly affects some of the outcomes of 
the analytical studies of preheating of bosons, which 
hardly go beyond the Hartree approximation. For this 
reason, the results presented in our work are obtained 
with numerical simulations on the lattice. More pre- 
cisely, the code "LATTICEEASY" l^, by G. Felder and 
I. Tkachev, has been used (details are given in section 
3). Full numerical calculations on the lattice are how- 
ever rather extensive. We have found that the neces- 
sary computing time is reduced in the conformal case, 
that is with the inflaton potential V {(p) = A0'*/4, and 
with a r.h. sneutrino mass which is negligible during the 
early stages of preheating. For this reason, in our com- 
putations we fixed Mn = 10^^ GeV, which is smaller 
than the Hubble parameter during inflation, but still high 
enough to require a nonthermal production of the sneu- 
trinos. The numerical results show a very efficient pro- 
duction of r.h. sneutrinos and inflaton quanta at preheat- 
ing/rescattering. Even for a coupling inflaton-sneutrino 
as small as ^ few x 10^^^ , the produced quanta come 
to dominate the energy density of the Universe already 
within about the first 5 e-folds after the end of inflation. 
In particular, the energy density stored in sneutrinos is 
typically found to be a fraction of order one of the total 
energy density, so that a sufhcient leptogenesis is easily 
achieved at their decay. 

R.h. sneutrinos are coupled to Higgs fields and left 
handed (l.h.) leptons through the superpotential term 
hN H L C W (responsible for the Yukawa interaction 
which provides a Dirac mass to the neutrinos). Thus, 
one may expect that quanta of the latter fields are ampli- 
fied by the rescattering of the r.h. sneutrinos produced 
at preheating. We study this possibility in section 4 , 
showing that indeed the amplification occurs for a wide 
range of values of the coupling h . Part of the analy- 
sis follows the detailed discussion on rescattering given 
in , where the numerical code \ ^7\ u sed here was also 
employed. However, the analysis of [|36[ is focused on the 
production of massless particles, while we show that the 
non vanishing mass of the sneutrinos can have some inter- 
esting consequences. More precisely, when the sneutrino 
quanta become non relativistic (let us denote by r} the 
time at which this happens) their rescattering effects be- 
come much less efficient. Thus, a strong amplification of 
the MSSM fields at rescattering can take place only if the 
coupling h is sufficiently large so that the amplification 
occurs before t) . As a consequence, for massive sneu- 
trinos and for small values of h , the number of MSSM 
quanta produced at rescattering is an increasing func- 
tion of h . However, the production is actually disfavored 
when the coupling h becomes too high. This is simply 
due to energy conservation, since the energy associated 
to the interaction term between the sneutrino and the 
MSSM fields cannot be higher than the energy initially 
present in the sneutrino distribution (equivalently, one 



can say that, for a too high coupling h , the non vanish- 
ing value of the sneutrinos gives a too high effective mass 
to the MSSM fields, which prevents them from being too 
strongly amplified). Posing quantitative bounds on the 
coupling h would require some better (analytical) under- 
standing of the details of rescattering than we presently 
have. However, the numerical results shown in section 4 
may give an idea of the expected orders of magnitude. 

An important remark is in order. When we speak 
about the amplification of MSSM fields coupled to the 
r.h. sneutrinos we have actually in mind amplification 
of _D— fiat directions (let us generally denote them by 
X). Indeed, D— terms provide a potential term of the 
form AV ~ for any scalar non flat direction Y . 

Since go is a gauge couphng {go = O (lO^^)), we expect 
such terms to prevent a strong amplification of Y , again 
from energy conservation arguments. Another important 
issue which emerges when gauge interactions are consid- 
ered is whether gauge fields themselves are amplified at 
rescattering. We believe that, at least in the model we 
are considering, also the amplification of gauge fields will 
be rather suppressed. The scalar distributions amplified 
at rescattering break much of the gauge symmetry of 
the model. This gives the corresponding gauge fields an 
effective mass in their dispersion relation (analogous to 
the thermal mass acquired by fields in a thermal bath) 
of the order ^ gQ{X'^) . As we extensively discuss 
in the paper, in the class of models we are considering 
the nonthermal distributions formed at rescattering are 
characterized by a typical momentum several orders of 
magnitude smaller than this mass scale. For this reason, 
one can expect that such heavy gauge fields cannot be 
strongly amplified. Jz?] In our opinion, an explicit check 
of these conjectures by means of numerical simulations 
could be of great interest, especially considering the great 
importance that gauge fields could have for the thermal- 
ization of the scalar distributions. 

To conclude, we discuss the production of gravitino 
quanta from the scalar distributions generated at rescat- 
tering. We already mentioned that in order to avoid a 
thermal overproduction of gravitinos an upper bound has 
to be set on the reheating temperature Tub of the ther- 
mal bath, Thh ^ few x 10^° GeV The requirement 
of a low reheating temperature can be seen as the demand 
that the infiaton decays sufficiently late, so that particles 
in the thermal bath have sufficiently low number densi- 
ties and energies when they form. If ~ 10^^ GeV at 
the end of infiation, and if the scale factor a is normalized 
to one at this time, the generation of the thermal bath 
cannot occur before a ~ 10^ . Gravitino overproduction 
is avoided by the fact that in the earlier times most of 
the energy density of the Universe is still stored in the 
coherent infiaton oscillations. On the contrary, we have 
already remarked that preheating/rescattering lead to a 
quick depletion of the zero mode in the first few e-folds 
after the end of infiation. |^ 

The question whether also the distributions formed at 
rescattering may lead to a gravitino problem is thus a 
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very natural one, and section 5 of the paper is devoted 
to some considerations on this regard. [£9| To provide 
at least a partial answer to this question, we distinguish 
the period during which rescattering is actually effective 
from the successive longer thermalization era. The com- 
putation of the amount of gravitinos produced during 
the earlier stages of rescattering appears as a very diffi- 
cult task. The numerical simulations valid in the case of 
bosonic fields indicate that a perturbative computation 
(with dominant 2 — > 2 scatterings taken into account) can 
hardly reproduce the numerical results, and that prob- 
ably N ^ 2 processes {N > 2) have also to be taken 
into account (we discuss this point in more details in sec- 
tion 4). It is expected that the same problem will arise 
also for the computations of the quanta of gravitinos pro- 
duced by the scalar distributions which are being forming 
at this stage. The end of rescattering/beginning of the 
thermalization period is instead characterized by a much 
slower evolution of the scalar distributions. In particu- 
lar, the total occupation number of all the scalar fields is 
(approximatively) conserved, which is interpreted |^ by 
the fact that 2^2 processes are now determining the 
evolution of their distributions. Motivated by this obser- 
vation, we assume that 2 — > 2 interactions are also the 
main source of production for gravitinos from this stage 
on. 

In the thermal case, the gravitino production is dom- 
inated by processes having a gravitationally suppressed 
vertex (from which the gravitino is emitted) and a second 
vertex characterized by a gauge interaction with one out- 
going gaugino. However, we believe that in the present 
context these interactions will be kinematically forbid- 
den, due to the high effective mass-squared that gaug- 
inos acquire from their interaction with the scalar dis- 
tributions (the argument follows the one already given 
for gauge fields). Once again we notice that the system 
is still effectively behaving as a condensate: the num- 
ber densities of the scalar distributions are set by the 
quantity {X^) , which is much higher than the typical 
momenta of the distributions themselves. This generates 
a high effective mass for all the particles "strongly" cou- 
pled to these scalar fields. A further comparison with 
the case of a thermal distribution may be useful: in the 
latter case both the typical momenta and the effective 
masses are set by the only energy scale present, namely 
the temperature of the system. As should be clear from 
the above discussion (see also [Q), the thermalization 
of the distributions produced at rescattering necessarily 
proceeds through particle fusion. Only after a sufficiently 
prolonged stage of thermalization, will the system be suf- 
ficiently close to thermodynamical equilibrium so to ren- 
der processes as the one discussed above kinematically 
allowed. 

In section 5 we show that if this class of processes is 
indeed kinematically suppressed, the production of grav- 
itinos from the distributions formed at rescattering is 
sufficiently small. However, we remark that this anal- 
ysis still leaves out the gravitino production which may 



have occurred at the earlier stages of the rescattering pe- 
riod. Whether this production may be sufficiently strong 
to overcome the limits from nucleosynthesis remains an 
open problem. 

Let us finally summarize the plan of the paper. In 
section 2.1 wc introduce our notation and briefly discuss 
some neutrino low energy phenomenology. In section 2.2 
we discuss leptogenesis with a thermal production of the 
r.h. (s)neutrinos. Leptogenesis with a nonthermal pro- 
duction of r.h. neutrinos is reviewed in the following 
subsection. The supersymmetric version of this model is 
presented in section 3 , where we study the nonthermal 
production of sneutrino quanta. Section 4 is devoted to 
the amplification of the MSSM _D— fiat directions due to 
the rescattering of the r.h. sneutrino quanta. The discus- 
sion on the gravitino production is presented in section 
5 , apart from a few technical details which can be found 
in the appendix. 



II. SEE-SAW PHENOMENOLOGY AND 
LEPTOGENESIS 

In this section, we introduce our notation for the SUSY 
see-saw and outline its low-energy implications. The aim 
is to make contact between realistic see-saw models, and 
the one generation toy models in which we will study 
the sneutrino production. We discuss the lepton asym- 
metry that can be produced in (s)neutrino decay, which 
implies a lower bound on the mass of the lightest r.h. 
(s)neutrino. Then, we briefly review different mecha- 
nisms for r.h. neutrino production, namely thermal and 
non thermal. The terms neutrino and sneutrino will be 
used interchangingly in discussing thermal production, 
which is similar for bosons and fermions. Concerning the 
nonthermal case, instead, differe nt re sults are obtained 
for the two species, and in section II C we review the ones 
for the neutrinos. Nonthermal production of sncutrinos 
is instead discussed in the next section. 



General Framework 



Let us consider the Minimal Supersymmetric Standard 
Model (MSSM) extended with three r.h. neutrino super- 
fields Ni (sometimes called the minimal supersymmetric 
see-saw model). The relevant couplings of the r.h. neu- 
trinos are given by the superpotential 



Wn 



hji Li 



1 



(1) 



where Li and H are the lepton and the Higgs doublets, 
respectively, and /i is a 3 x 3 complex Yukawa matrix. We 
will neglect the phases in our analysis of N production, 
because CP violation is not required for this process. 
We work in the r.h. neutrino mass basis, where the mass 
matrix M is diagonal, and we disregard the possibility of 
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Mechanism 


N Yukawa h 


N mass 


(j>-N coupling 


Thermal 


IQ-^oV < mi < 10^ 




lO'-'GeV < Afi < Trw 


irrelevant 


Affleck-Dine 


IQ-^eV < < 10 






f < //infl 

i (Mr) < 


Pert.0 1 
decay J 


Tlv < H{Ti) 




f < 771^/2 1 
iWi > / J 


Bi?(0 ^ A.A,) ~ 1 

JD Tl,[(p — ' iVj iVj ) ~ 1 


A'' preheating 
eq. (@ 


Vlv < H{n) 




Mi > 10" GeV 


> 0.03 


A*' preh./resc. 
eq. (PI) 


Vlv < H{n) 




M, ^ g, 10^^ GeV 





TABLE I: Summary ol parameters for which leptogenesis could work, for different r.h. (s)neutrino production mechanisms. 
In (s) neutrino production mechanisms which do not require the Yukawa coupling (non-thermal mechanisms), the constraint 
on the Yukawa matrix is that lepton number violating interactions in the thermal soup be out of equilibrium after the r.h. 
(s)neutrinos decay at n. This also implies Mi > T{Ti), and possibly additional constraints on L violating processes mediated 
by Nj,j 7^ i. Recall fhi parametrises the Ni decay rate, and is defined after eqn. ( |l2[) . The Affleck-Dine mechanism proceeds 
through generation of large expectation values either for a small or a tachyonic pH effective mass of sneutrinos M'^^ during 
infiation. The asymmetry made by the perturbative decay of the infiaton can be generated by the on-shell r.h. (s)neutrinos [ p6[ , 
which subsequently decay, or by the decay via off-shell r.h. (s)neutrinos (A*) to Higgses and leptons |^]. The properties of 
nonperturbative (s)neutrino production analysed in the present paper is summarized in the last two lines (A is the infiaton self- 
coupling). Other scenarios for nonperturbative production after the end of inflation can be envisaged, with model-dependent 
results. 



nearly degenerate r.h. (s)neutrinos ||Tl|, gs], (i.e. we 
assume that the difference of neutrino masses is of order 
their mass). 

The lepton asymmetry produced in the decay of Ni 
can be written 



Yr = 



Nr. - Nr 



N, 



Ni 



(2) 



where N^. is the total number density of the zth heavy 
(s)neutrino species prior to its decay, s is the entropy 
density at decay , Hi parametrises washout effects due 
to subsequent lepton number violating interactions, and 
Ci arises from the CP violation of the Ni decay. It is 
given by [^l 



r{N, 



LH)-T N,, 



r(iv. 



8^(Mt)^ 



LH)+T(N, 



^Im (M^) 




(3) 



where f{x) = y/x[2/{x— 1) +ln(l/.T + 1)] for hierarchical 
r.h. neutrino masses. 

We suppose for the moment that some number density 
of Ni is produced in the early Universe, and concentrate 
on how large an asymmetry can be generated. The asym- 
metry ei is determined by the masses and couplings of the 
r.h. (s)neutrinos, which are given in eqn. (|l|). However, 
it can be related to, and therefore constrained by, low 
energy observables. 



The CP asymmetry produced in the decay of a r.h. 
(s)neutrino can conveniently be parameterized as 



it 



3 mg 



Sep 



M,, 



IQiOGeV 



VO.OSeV/ 



'CP- 



(4) 



By using eqs. (|^) and (||), it is possible to show ||Tj^ 
that for the case ei, 5cp satisfies the upper bound~p2| 




\5cv \ < 1 



(5) 



By combining the two last expressions, one finds an upper 
bound on the parameter ei which scales linearly with 
the r.h. (s)neutrino mass Mi . We will shortly see that 
this implies a lower bound on Mi for leptogenesis to be 
viable. The mass in equation (U) denotes the mass 
of the heaviest left-handed neutrino. The light neutrino 
mass matrix is obtained by integrating out the heavy r.h. 
neutrinos to give the see-saw formula 



0\2 



(6) 



We will assume that the light neutrino masses rm are 
hierarchical, so TO3 ~ y^Amf^ 

If h is written in the charged lepton mass eigenstate 
basis (neutrino flavour basis), then rrii, is diagonalised 
by the MNS matrix U [Q, which can be written U = 
V ■ diag(e-^'^/2^ e-^-^'/^, 1), with 
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V 



C13C12 C13S12 Sise""' 

-C23S12 - S23Sl3Cl2e*'' C23C12 - S23Sl3Sl2e"' S23C13 
S23S12 - C23Sl3Ci2e*'^ -S23C12 - C23Sl3Sl2e"' C23C13 



In this matrix, C12 = cos0i2, and so on. Atmospheric 
I, in and solar [||, || data imply that 6*12 and 6*23 
are large, approaching 7r/4. ^13 is constrained to be 
< 0.1 by the CHOOZ experiment |44|. In a supersym- 
metric scenario, there is additional information about h 
and M available in the slepton mass matrix. The neu- 
trino Yukawa h'^h appear in the renormalization group 
equations for the soft slepton masses, and thereby in- 
duce flavour violating slepton mass terms [m|]y . 
In a simple-minded leading log approximation, these off- 
diagonal mass matrix elements are 



(3m^ + Al) 



MgutJ 



(7) 



where hi are the eigenvalues of h, mg and Aq are soft 
parameters at the GUT scale, and we introduce a new 
matrix Vl which diagonalises /i^/i in the charged lep- 
ton mass eigenstate basis {VLh''hVl = diagonal). The 
branching ratio for £j can be roughly estimated as 

i3: 



■ tan^ f3 



(8) 



where m| is the slepton mass scale. The experimen- 
tal bound BR{^i 67) < 1.2 x 10"" iQ implies 
[jTilJ^e ^ lO'^rhj^, for — 100 GeV. This constrains 
the angles in Vl for given neutrino Yukawas hi. It can 
be shown that M and h have the same number of pa- 
rameters as the weak scale neutrino and slepton mass 



matrices. Furthermore, in a SUSY scenario with uni- 
versal soft masses at the GUT scale, M and h can be 
parametrised with ifij^ and nii, [ |7| . The r.h. neutrino 
masses and Yukawa couplings can be therefore recon- 
structed {in principle, but not in practise |47 ) from the 
weak scale neutrino and sneutrino mass matrices, so that 
Ei can be expressed in terms of weak scale variables. An 
analytic approximation for ei can be found in Esl Eol : 



ihl 



Im 



(9) 



where rrii are the light neutrino masses, hi is the smallest 
eigenvalue of h, and W = VlU is the rotation from the 
basis where the i^l masses are diagonal to the basis where 
h^h is diagonal, hi is in practise unmeasurable; however, 
if h has a hierarchy similar to the up Yukawa matrix 
hu, then h^ ~ 10~®, and ei will only be large enough 
if there is some enhancement from the imaginary part. 
There are two simple limits for the matrix W, which are 
motivated by model building. The first is Vl — 1, and 
corresponds to an almost diagonal slepton mass matrix 
(in the charged lepton mass eigenstate basis) . This means 
that the large mixing observed in the MNS matrix U 
must come from the r.h. sector [Q. The second option 
is W ~ 1, so Vl — C^^- This would arise if the large i^l 
mixing is induced in the l.h. sector ||5^. In the Vl — 1 
case, eqn. (P) gives [Esl: 



ei 



Im ■ 



TO? c?,c 



13''12 



3hl 
'SttD 



"^2 C?3S?2 



mi si, e^^' 



TOi 0^30^2 e*"^ -I- TO2 cfgsfj e*"^' -I- TO3 sfg 



b') - sin(0 - ^') 
m2 



p2i5 



(10) 

(11) 



where D — m\ C13C22+TO2 Ci3Si2+w| 5^3, and in the sec- 
ond equation, the solar and atmospheric angles have been 
taken to be 7r/4. If we estimate the phases to be 0(1), 
hi ^ the up Yukawa, and TO3/TOI ~ Ato^^^/Ato^^;, this 
gives e < 10~^(si3/.l)^, where we have scaled the unmea- 
sured angle 6*13 by its upper bound. This is barely 
large enough for thermal leptogenesis. However, we re- 



mind that hi is unknown and it can well be hi > 10" . 
The second case, where ~ 1, can arise if M and h^h 
are almost simultaneously diagonalisable |Q. For small 
angles in the approximation for e can be extracted 
from (^0|), replacing the angles of the MNS matrix by 
the angles of W, and setting the cosines — > 1. When 
W ^ 1, then Vl so it is the MNS angles that ap- 
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pear in equation (0), and BR{^i e-f) < 1.2 x 10"" @ 
implies an upper bound on the CHOOZ angle 6*13 < .02 
(for mL = 100 GeV, ^13 = 1)[||. To conclude, we briefly 
comment on the parameter k^. After the asymmetry is 
generated in the out-of-equilibrium decay of the r.h. neu- 
trino, lepton number violating interactions which could 
wash out the asymmetry must be out of equilibrium. 
This is a fairly straightforward requirement when con- 
sidering the decay of the lightest r.h. neutrino A^i [ p^ ; 
it is more complicated in the case of iV2.3 decaying at 
T > Ml . The fraction of the asymmetry which sur- 
vives these interactions is < 1. 



Ml > (10^ - 10^°) GeV, although it is fair to say that 
higher values are required if the bound ^ is not satu- 
rated. If A^i is generated from the thermal bath, a re- 
heating temperature greater than Mi is required. In the 
most favorable case, the required value is only marginally 
compatible with the bound imposed by gravitino over- 
production from the thermal bath p^ . To overcome 
the potential conflict between the gravitino bound and 
the requirement of a reheating temperature high enough 
for leptogenesis, the possibility of producing right-handed 
neutrinos non thermally has been envisaged . The fol- 
lowing considerations will be focused on this framework. 



B. Thermal A'^; production 

We now consider the case where the lightest r.h. 
(s)neutrino A^i is thermally produced after Trh [^sj. 
With hierarchical r.h. neutrino masses, one can typi- 
cally assume the lepton asymmetry to be produced by 
the decay of the lightest r.h. (s)neutrino A^i . As we shall 
see, this is a self-consistent assumption, because M2 and 
M3 will turn out to be larger than Tr//. To generate a 
lepton asymmetry, the decay of the Ni should proceed 
out of equilibrium. More quantitatively, the ratio of the 
thermal average of the Ni decay rate and of the Hubble 
parameter at the temperature T ~ Mi , 



2fj\T=M. - 2 X 10-3 eV 



(12) 



should be less than unity to have an unambiguously 
out-of-equilibrium decay. The parameter ifii is defined 
as Ti{H)'^ /Mf , where (H) is the Higgs vev. However, 
mi cannot be taken too small if A^i is produced ther- 
mally ||l2| . Indeed the quantity rhi controls the strength 
of the interactions of the A^i with MSSM degrees of free- 
dom, and an efficient thermal production via Yukawa in- 
teractions typically requires rhi > 10~^ eV. To account 
for both these effects, Tn^ ~ H{T = Mi) should be 
taken, so the decay is only barely out of equilibrium, 
and the final lepton asymmetry has to be computed by 
integrating the full set of relevant Boltzmann equations 
pO| , [Tlj l . These computations show that a significant por- 
tion of the lepton asymmetry is erased by lepton-number 
violating processes, and that only a fraction k < .1 or 
less typically survives. Starting from A^i in thermal equi- 
librium at T > All , and collecting all the above infor- 
mations, the final baryon asymmetry can be estimated 
to be 



Yn 



10" 



-10 



Ml 



101" GeV 



)( 



0.05 eV 



CP 



(13) 



This expression has to be compared with the baryon 
asymmetry required by Big Bang Nucleosynthesis Yb = 
Nb/s ~ (1.7 - 8) X 10"" i, 11. As we have antici- 
pated, we see that thermal leptogenesis can be a viable 
mechanism if the mass of the r.h. neutrino Ni is suf- 
ficiently high. From eqn. ( p^ ) we find the lower bound 



C. Non thermal production of right handed 
neutrinos 

Many alternatives to thermal (s)neutrino production 
have been considered in the literature. We concentrate 
here on mechanisms that involve a direct coupling of N 
to the inflaton, although in the next section we will com- 
ment on differences and similarities with the Affleck-Dine 
mechanism. The strength of this interaction, relative to 
the inflaton coupling to other degrees of freedom, is a free 
parameter; for appropriate values, a lepton asymmetry 
of the correct magnitude can be produced. The number 
density of r.h. (s)neutrinos will also depend on the evo- 
lution of the inflaton between the end of inflation and 
reheating. If the inflaton decays perturbatively, right- 
handed neutrinos with masses less than half the inflaton 
mass could be produced in the decay |Q. For heavier 
r.h. neutrinos, one can also envisage the possibility that 
a sufficient leptogenesis is generated in processes in which 
they mediate a perturbative inflaton decay \ ^7\ . In both 
cases, the final lepton asymmetry will be proportional 
to the branching ratio of the inflaton into (either on- or 
off-shell) neutrinos. A branching ratio of order one is 
typically required. Right-handed neutrinos with masses 
greater than that of the inflaton can be produced at pre- 
heating, if their interaction with the inflaton is strong 
enough. The production of heavy fermions (sneutrinos 
are discussed in the next section) in an expanding Uni- 
verse was first discussed in ref. [g^ (fermionic production 
in the conformal case was first studied in Jssf), where a 
direct Yukawa coupling to the inflaton cj) was considered, 
and the simplest chaotic inflationary scenario with a mas- 
sive inflaton, V {(j)) = m|(/)^/2 , ~ 10^^ GeV, was 
assumed. The relevant part of the lagrangian is 



Cn,4> = N (M + g^) N , 



(14) 



where N is any one of the r.h. neutrinos. We assume that 
only one r.h. neutrino generation plays an important role 
in the generation of a lepton asymmetry, and therefore we 
drop the r.h. neutrino generation index for the remain- 
der of this section. The generalization of the following 
analysis to three generations is straightforward, at least 
as long as the r.h. neutrino-inflaton coupling matrix g 
is diagonal in the r.h. neutrino mass basis (otherwise, 
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the formahsm of |Q should be used). After the end of 
inflation, the inflaton condensate (/> oscillates about the 
minimum of its potential with amplitude of a fraction 
of the Planck mass Mp ~ 1.22 • 10^^ GeV . The total 
effective mass of the fermion M + g (f){t) varies non adia- 
batically in time, and this leads to a (non perturbative) 
production of quanta of N. In particular, fermion pro- 
duction at preheating occurs whenever the total effective 
mass crosses zero. As a consequence, fermions with a 
mass up to 



10 



10 



1/2 



lO^^GeV, q . 



9^ 
4m?. 



35^ loio 



_ (15) 
can be produced |3J] , irrespective of the value of the re- 
heating temperature of the thermal bath which is formed 
at later times. The abundance of neutrinos produced at 
preheating has been computed analytically ||35|] , and it is 



most conveniently given in terms of the ratio 



P0 



1 



1.4 X 



10-" g 



1010 GeV 



M- 



1/2 



10 



log 1.7 X 10-^ 



9 



1/2 



3/2 



Ai'lO. 

(16) 

where we have defined Mio — M / (10^° GeV) . The above 
formula is valid as long as the backrcaction of the pro- 
duced neutrinos on inflaton dynamics is negligible, as it 
turns out to be the case as long as g ^ 10* |3^. For 
larger values of g, the effectiveness of preheating increases 
(by a factor up to about 1.5), and the above equation 
gives just a lower bound on A^^r. In what follows, we 
will conservatively assume q < 10* . For a massive 
inflaton, the ratio (|6|) is constant until the inflaton con- 
densate decays. If neutrinos decay before reheating has 
completed, the resulting baryon asymmetry reads 



Yn 



23^^-^ 
4 -10-" 



23 



-eA^, 



3 T, 



N 



RH 



IQ8 10 



1/2 



4 



109 GeV 0.05 eV 
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Log 1.66 • 10 



3 1 



1/2 



Al- 



io 



3/2 



(17) 



The ratio of Nn to the entropy is constant after reheat- 
ing, unless it decreased due to some subsequent entropy 
production For a given value of the parameter g, the 
baryon asymmetry (18) is maximized when the mass M 
has its largest possible value (16). For Trh — 10^ GeV, 
m^^ = 0.05 eV and 6cp — 1, imposing the condition 
(16) on eqn. (18) it is possible to see that the genera- 
tion of the observed baryon asymmetry requires g larger 
than about 0.03. In particular, if we assume q ~ 10*, 
then the observed baryon asymmetry can be obtained 
for a mass of the right-handed neutrino of the order of 
10^^ — 10^^ GeV. Production of r.h. neutrinos at pre- 
heating can generate a large enough lepton asymmetry 
from neutrinos with Mjy > Trh , so constraints on Tpji 
do not translate into bounds on M^. We notice how- 
ever that the large inflaton-neutrino coupling required 
{g > 0.03) can in principle modify through quantum ef- 
fects the small mass or self coupling parameters charac- 
terizing the inflaton potential. This constitutes a further 
motivation for considering supersymmetric models, as we 
do in the remaining of this work. We will see that the 
production of the scalar partners of the r.h. neutrinos 
can significantly affect some of the above considerations. 



III. SUPERSYMMETRIC SEE-SAW AND 
NONTHERMAL PRODUCTION OF RIGHT 
HANDED SNEUTRINOS 



As we have seen in the previous section, preheating 
can have very important consequences for leptogenesis 
through the production of r.h. neutrinos [^^. In super- 
symmetric extensions of the see-saw model, the produc- 
tion of the supersymmetric partners of the neutrinos is 
even more important. Due to supersymmetry, the in- 
flaton couples with the same strength both to the r.h. 
neutrinos and to the sneutrinos, so that if the former are 
produced at preheating this will also occur for the latter. 
However, while production of fermions is limited by Pauli 
blocking, the production of scalar particles at preheating 
is characterized by very large occupation numbers. As a 
consequence, the production of r.h. sneutrinos can be ex- 
pected to be more significant than the one of neutrinos, as 
the numerical results presented below confirm. Produc- 
tion of particles at preheating gives very model depen- 
dent results; nevertheless, some general features can be 
outlined, and the whole process can be roughly divided 
into three separate stages. The first of them is charac- 
terized by a very quick amplification of the fields directly 
coupled to the inflaton (and of the inflaton fleld itself, 
in the case of a sufficiently strong self-interaction) to ex- 
ponentially large occupation numbers ]26t . Very rapidly, 
the system reaches a stage in which the backrcaction of 
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the produced quanta, customarily denoted as rescatter- 
ing |3l|, plays a dominant role. In the case of parametric 
resonance, the scatterings of the quanta against the zero 
mode of the inflaton destroy the coherence of the oscilla- 
tions, thus ending the resonant production characterizing 
the early stage of preheating P6[ . An equally important 
backreaction effect is the amplification of all the other 
fields to which the produced quanta are coupled. This is 
a very turbulent process, dominated by the nonlinear ef- 
fects caused by the very high occupation numbers of the 
fields involved. As a result, all these mutually interacting 
fields are left with highly excited spectra far from thermal 
equilibrium. The latter is actually achieved on a much 
longer timescale, through an adiabatic (slow) evolution of 
the spectra, which characterizes the third and final stage 
of the reheating process. ||8^ The first stage of preheating 
is well understood. Particle production is computed in 
a semi-classical approximation (for a rather general for- 
malism in the case of several coupled fields see ) , and 
analytical solutions have been obtained in a broad class of 
models ||2^, |6[ ^ |5^, |6^ . Analytical approximations 
break down when nonlinear processes become dominant. 
However, the high occupation numbers of the scalar fields 
involved allow a classical study of the system. Indeed, in 
the limit of high occupation numbers quantum uncertain- 
ties become negligible, and quantum probabilities show 
a classical (deterministic) evolution ||6^. The latter can 
be better computed by means of lattice simulations in 
position space |3^ , where all the effects of back- 

reaction and rescattering are (automatically) taken into 
account. A detailed discussion of rescattering and of the 
approach to thermal equilibrium has been given in 
where the code "LATTICEEASY" by G. Felder and 
I. Tkachev, has been used. The numerical results pre- 
sented in this paper are also obtained with this code. 
For numerical convenience, we consider a chaotic infla- 
tionary scenario with a quartic potential for the inflaton. 
More specifically, we focus on the superpotential [p9[ 



W^($,7V) 



1 



(18) 



The second term of W reproduces the lagrangian (|lj) 
for the r.h. neutrinos. We denote the scalar components 
of the inflaton and of the r.h. neutrinos multiplets with 
(j) and N , respectively. To simplify the numerical com- 
putations, the imaginary components of the scalar fields 
will be neglected. Therefore, after canonical normaliza- 
tion, (j) -^6/^/2, N — > we consider the scalar 
potential |]90| 

Scalar = ^ 0* + ^ (ff + Mf ■ (19) 

The size of the temperature fluctuations of the Cosmic 
Microwave Background sets A ~ 9 x 10"^"* , while the 
neutrino mass M as well as the coupling g to the inflaton 
are model dependent parameters. The case M = is 
analyzed in detail in [ p6| . In figure |l| we show the time 
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FIG. 1: Time evolution of the comoving number density 
and of the effective masses of the inflaton cj> and of the r.h. 
sneutrino A*' . See table ^ for notation. 



evolution of the comoving number density and of the co- 
moving effective mass of the two scalars and iV , for 
the particular choice of the parameters M = lO^^GeV, 
g = g^/A = 200 . The effective mass is defined as 



= a2 



a 
a 



(20) 



where a is the scale factor of the Universe, normalized to 
one at the end of inflation, prime indicates derivative with 
respect to the conformal time 77 , while (. . .) denotes aver- 
age over the sites of the lattice. The term a" /a appears in 
eqn. (^ ) because we are considering minimally (rather 
than conformally) coupled scalars, and it vanishes in a 
radiation-dominated background. The comoving num- 
ber density is defined as the integral over momentum of 
the "occupation number" 



fc2 



1 



cff ^ 



(21) 



where fk denotes the Fourier transform (to be evalu- 
ated on the lattice) of the rescaled field a <j) . By defi- 
nition IstII , the quanta stored in the oscillating inflaton 
condensate do not contribute to Nc^^ in figure |l]. The 
three quantities rricff , , and rj are all shown in units of 
V\(j)o o± 1.25 • 10^2 Qgv, with 00 - 0.342 Afp denoting 
the value of (f> at the end of inflation, to the appropri- 
ate power. All the numerical results presented in this 
work are obtained ^with a two dimensional lattice of size 

i = 20 (yX<j)o) and with N = 1024^ sites (see fs^ 

for details). Figure]^ exhibits the features that we have 
outlined at the beginning of this section, namely a quick 
stage of exponential growth of the occupation numbers 
followed by a period in which the occupation numbers 
are nearly constant. During the first stages of the pro- 
cess, the results presented reproduce very well the ones 
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inflaton vev at the end of inflation 
inflaton self-coupling (eqn. (Q) 
neutrino- inflaton coupling (eqn. (|l^)) 
neutrino Yukawa coupling (eqn. (p6[)) 
X number density (eqn. (pl|)) 
X "mass" (eqn. (poh) 



(Po ^ Mp/3 
A ~ 9 X 10-" 
g = //A 



h = h^/\ (* 



Nc,x ~ [comov. num. den.]/[%/A<^o]^ 



- ^) /[^/A<^o]^ 



TABLE II: Translation table between quantities in plots and superpotential parameters (we give the eqn. where the parameter 
is defined). Recall that -q is the conformal time coordinate, 77 = at the end of inflation, and subsequently the scale factor is 
a{-n)/a{0)~r,/2 + l. 



obtained in for M = , since the "bare" mass of 
the r.h. neutrinos is initially negligible. However, the 
presence of a non vanishing bare mass affects the subse- 
quent evolution of the system. Indeed, when the value 
of the fluctuations in the sncutrino field becomes com- 
parable with the amplitude of inflaton oscillations, the 
second term in eqn. ( p^ ) shifts the minimum of the effec- 
tive potential of the inflaton, giving it an effective mass 
that is roughly constant in comoving units. A stronger 
effect is related to the fact that the sneutrino itself is 
massive. In rescaled units and for the present choice of 
the parameters, we have 




0.08 a 




(22) 



Numerical results show that (ip) ~ <po/a (as one could 
also see by inspecting the potential (|l3) in Hartree ap- 
proximation), so that the part in the above equation 
that depends on (f> remains of order one. It follows 
that the r.h. neutrino mass M dominates the effective 
mass ( p^ for > 77 ~ 350 , as clearly indicated by the 
growth of mcff.N visible in figure ^ for r/ > ?) . Produc- 
tion of sncutrinos at preheating in the present model is 
strictly related to the production of fermions we have 
analysed in section 2.3 . In particular, production oc- 
curs whenever the effective neutrino mass ( p^ crosses 
zero. Numerical results show that preheating is termi- 
nated by rescattering effects when the scale factor a is 
of the order of Oresc — 100. As a consequence, sncutri- 
nos with a bare mass up to g(j)o/arcsc — g ■ 10^^ GeV 
will be efficiently produced at preheating and will con- 
stitute a sizable fraction of the background energy [ pl| . 
Numerical results show that after the onset of rescatter- 
ing, the energy density gets roughly equiparated between 
the quanta of the two species. As a consequence, in gen- 
eral large couplings correspond to a large interaction en- 
ergy, and therefore to a smaller number density during 
the rescattering/thermalization stage In partic- 

ular, for this reason a large quartic self-coupling \N\* 
for the sneutrino would prevent it from getting large oc- 
cupation numbers, since energy conservation would im- 
pose (iV^) cx \/X/g^- Numerical results also indicate 
that soon after the beginning of rescattering most of the 
r.h. neutrino quanta have a momentum of the order 
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FIG. 2: Time evolution of the derivative of the scale factor 
with respect to conformal time. A constant value indicates 
radiation domination, while a linear growth indicates matter 
domination. 



/c* 15 \/X(j)o ■ ||9|] Thus, most of the r.h. neutrinos be- 
come non relativistic at a time not much greater than ?). 
From this time on, the energy density of the system red- 
shifts as the energy density of matter. The transition be- 
tween the two stages of matter and radiation domination 
is clearly visible in figure H, where we show the time evo- 
lution of the derivative of the scale factor with respect to 
conformal time. As long as the neutrino mass is negligi- 
ble, the energy density of the system redshifts as the one 
of radiation ^ , and the evolution of the scale factor 

is very well approximated by a ~ (1 + 1)^^^ — ri/2 + 1 , 
where we have set i = 77 = at the end of inflation. 
Therefore, during the initial stage of radiation domina- 
tion, a' is constant. In the following matter dominated 
stage a oi rf , and a' grows linearly with time. To esti- 
mate the baryon asymmetry produced from the decay of 
the r.h. sneutrinos, we need to know the fraction of the 
entropy of the Universe that is generated in the decay. 



The baryon asymmetry will be |13 



^ ~ '23 SN stot ~ 23 I ^4 M ' 
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where s^r cx g^,{Tpf)T^ is the entropy produced in N de- 
cay, and Stot is the total entropy of the Universe. If the 
r.h. sneutrinos dominate the Universe when they decay, 
then T/v is the temperature to which the Universe is re- 
heated by the decay of the sneutrinos A'', and stot ^ sn- 
This condition is satisfied if the inflaton mainly decays 
only into one r.h. sncutrino family (as it is clearly the 
case in the one generation model we have studied numer- 
ically). However, if the inflaton couples to other scalars 
(also in particular to the other generations of sneutrinos) , 
these could produce additional entropy. At variance with 
the case of leptogenesis induced by the decay of right- 
handed neutrinos, analysed in section 2.3, the baryon 
asymmetry ( |2^ ) does not depend on the r.h. (s)neutrino 
mass, that must only satisfy M > Tjsf in order to prevent 
thermal regeneration of the r.h. sneutrinos after their de- 
cay. This is due the fact that, thanks to Bose statistics, 
r.h. sneutrinos can get large occupation numbers at pre- 
heating (whereas Pauli blocking makes fermion produc- 
tion less efficient), and they can easily represent a sub- 
stantial fraction of the energy in the Universe. |93 The 
generalization to the more realistic case of three neutrino 
families coupled to the inflaton is obtained by consider- 
ing the following supcrpotential, in the mass eigenstate 
basis for the r.h. neutrinos, 

W=^^^ + hj, U ■ Hu + i {M,6jk + V2g,fe$)7Vfe TV, . 

(24) 

This gives a potential for the real components of the 
scalars 



V = 



(25) 

where dots include the terms involving the Yukawa h, 
which are not relevant for nonthermal Ni generation. We 
suppose for simplicity that the neutrino-inflaton coupling 
gij is diagonal. Energy considerations after rescatter- 
ing ^ lead to the expectation that the sneutrino 
family that is most strongly coupled to the inflaton is 
also the one that will have the smallest number density 
(clearly, on the assumption that all the sneutrinos are 
sufficiently coupled to be amplified). This is opposite to 
the scenario in which sneutrinos are produced by the per- 
turbative decay of 0, where iVjVi is proportional to the 
inflaton branching ratio to Ni. However, the presence of 
the r.h. sneutrino bare masses, as well as the existence 
of couplings to the Standard Model degrees of freedom 
(whose effects will be analysed in detail in the next sec- 
tions), can strongly affect these conclusions. Although 
the resulting baryon asymmetr y f^ ) has the same ex- 
pression as the one reported in |13|, the mechanism that 
led to a sneutrino dominated Universe is different from 
the generation of large expectation values considered in 
ref. [p3| or in the Affleck-Dine [Bsl mechanism. Indeed, 



the latter is effective if during inflation the sneutrino (or, 
more generally, the amplified flat direction) has a mass 
much smaller than the Hubble rate. This requires (be- 
sides a sufficiently small bare mass M ^ H) that the am- 
plified field does not get a large effective mass through 
its coupling to the inflaton. It is important to remark 
that the mechanism we are discussing can provide a suf- 
ficient leptogenesis even if the coupling g between the 
inflaton and the r.h. neutrino multiplets is much smaller 
than the one needed in non supersymmetric models, i.e. 
with only the production of the neutrinos taken into ac- 
count, see eqn. (17). Anyhow, even couplings as small 
as g^ ~ 10 A prevent the formation of a large conden- 
sate during inflation. Therefore, the two mechanisms 
can lead to large occupation numbers in complementary 
regions of the parameter space. Notice that the above 
discussion applies to every effective mass term that can 
arise in the potential for a (quasi) flat direction. In par- 
ticular, it could be interesting to consider the effective 
mass of the order of the Hubble parameter that is gen- 
erally induced by supergravity corrections , although 
in this case ampliflcation effects may be weakened by 
the quick redshift characterizing the nonrenormalizable 
interactions. If supergravity corrections induce a tachy- 
onic mass m^g ~ —H"^, a large expectation value will 
be generated during inflation |24| , and the dynamics of 
preheating will turn out to be rather different from the 
one considered in the present section. This however re- 
quires a suitable nonminimal Kahler potential, and we 
will not discuss this possibility in this work. Finally, it 
is worth stressing that both the leptogenesis scenario de- 
scribed in this section and the one considered in [p^ , 
although they are related to Affleck-Dine leptogenesis, 
are somehow different from it for what concerns the ful- 
fillment of the Sakharov CP-violation condition . In 
the Affleck-Dine scenario, the latter is achieved by the 
motion of the Affleck-Dine condensate (that requires co- 
herence over many Hubble lengths), while in the mech- 
anism we are analysing, CP-violating sneutrino decays 
are crucial in the generation of an asymmetry. 



IV. PRODUCTION OF LIGHT DEGREES OF 
FREEDOM AT RESCATTERING 

The description presented above is further modifled 
by the effects of the coupling of the r.h. neutrino mul- 
tiplets to the l.h. leptons and Higgses, coming from the 
supcrpotential (0). Considering for simplicity only one 
generation, the supcrpotential ( |l8| ) will be supplemented 
by 



AW = hNL ■ H 



(26) 



The resulting scalar potential contains several interaction 
vertices coming from F-terms. In addition, there are 
quartic contributions from the D-terms. The presence 
of the latter in the potential plays an important role for 
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our discussion, since rescattering mostly affects the D- 
flat directions. To see this, consider the case in which 
the left-handed selectron and the charged scalar Higgs 
vanish. The D-terms for the neutral Higgs and sneutrino 
then take the familiar form of the tree level MSSM Higgs 
potential, with replaced by 



9SV{2) 



9l 



(27) 



The directions that are not D-flat (i.e. the ones for which 
\vl \ 7^ |^^"| in the present example) are characterized by 
a large (gauge) quartic coupling in the scalar potential. 
Due to this coupling, they cannot be significantly ampli- 
fied by rescattering effects. On the contrary, D-flat direc- 
tions have only quartic couplings coming from F-terms as 
\dW ldN\^ , whose strength \t? will be typically taken ^ 1 
in all the cases considered below. Indeed these quartic in- 
teractions among the D-flat directions will be neglected 
in our computation, since they can be relevant only dur- 
ing the thermalization stage, when the variances of these 
fields have grown to be sufficiently large. The most im- 
portant F-terms arising from the total superpotential are 
of the form - \? \Nf (|£>lP -f |i?"P) . If we denote by 
X the D-flat combination Pl| = l^*^!; then the relevant 
coupling of X to the sneutrino field will be simply given 
by ft.^ A^^ X^ (as in the previous section, we consider for 
numerical convenience only real directions). Besides the 
quartic term cx I? X* , we will also neglect the interaction 
term oc (gc/) + M) N X'^ :^ M N X'^ , which is responsible 
for the late time decay of the r.h. sneutrinos (that is, the 
supersymmetric counterpart of the vertex which gives the 
decay of the r.h. neutrinos into Higgses and leptons). We 
thus consider the simplified model characterized by the 
scalar potential 



V (0, N,X)^-^^ + -{M + g 4>f N^ + -e X^ 



(28) 

Neglecting the imaginary directions of the scalar fields, 
as well as many of the interaction terms, allows a con- 
siderable reduction of the computing time needed for the 
numerical simulations (this is particularly welcome for 
the extensive computation that we discuss in the next 
section). The above discussion leads us however to be- 
lieve that the simplified model should well describe the 
main features of the preheating and rescattering process 
for the supersymmctrized see-saw model with a nonper- 
turbative production of the r.h. neutrinos. In figure ^ we 
show the time evolution of the comoving number den- 
sity of the quanta of X . As in the previous section, we 
have fixed M = 10" GeV, 5 = 200 , while different val- 
ues of the parameter h = jX are shown. Even if in 
the simplifled model ( |2^ ) the X field is not directly cou- 
pled to the inflaton, we see that (for suitable values of 
the coupling K) it can be highly amplified by the rescat- 
tering of r.h. sneutrino quanta. Figure || shows that the 
growth of number of X particles is roughly exponential in 
time. When the effective sneutrino mass is varying non- 
adiabatically in time and is not negligible with respect to 
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FIG. 3: Time evolution of the comoving number density of 
the light quanta X . See table ^ for notation. 



sneutrino typical momenta, the production of X particles 
cannot be analysed in terms of scatterings of sneutrinos. 
However, after the end of the parametric resonance pe- 
riod and the onset of rescattering, one can expect that 
a particle-like picture can give some description of the 
behavior of the system js^. In this case, if the domi- 
nant contribution to X production process were given by 
the 2 — > 2 scattering XX XX, the rate of growth 
of Nx should be proportional to h"^. The lattice results 
appear to present a milder dependence on h, suggesting 
that the NN XX scatterings alone cannot account 
for the production of X states. In a naive perturbative 
analysis, the contribution, from (m -I- 2) x iV XX 
processes, to the rate of growth of Nx, is proportional 
to the m X N ^ XX rate times a factor roughly given 
by Nf; / {Anp^) , where p - l5VX(f>Q/ a is the typical 
momentum exchanged. Due to the high density of sneu- 
trinos, the expansion parameter /i^ N"^ / (Airp^) is of order 
unity for the values of h we are considering. Therefore, 
strongly turbulent processes with many incoming sneu- 
trinos can contribute significantly to the rate of growth of 
Nx ■ This is confirmed by the fact that the total number 
of particles decreases during the stage of generation of 
the X states, thus showing that particle fusion processes 
are dominant at this stage . The main features shown 
in figure ^ are shared by the other evolutions with differ- 
ent g that we will consider in the next section, and they 
can be understood at least qualitatively. As could have 
been easily guessed, the timescale for the growth of Nc^x 
is a decreasing function of h . We also notice from the 
figure that the amplification of X becomes less efficient 
as the quanta of N become non relativistic, at 77 ~ 300 . 
This can be seen explicitly in figure ^, where we show 
the effect of the r.h. neutrino mass term on the growth 
of the comoving occupation numbers of the and X 
fields. If the two fields are both massless, the rescatter- 
ing of the quanta of N lifts the X to (practically) the 
same amplification |p6| . We observe that the situation is 
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FIG. 4: Comparison of the growth of the occupation numbers 
of A*' and X with and without the r.h. neutrino mass term. 
Notice that the amplification of X considerably weakens when 
the r.h. neutrino quanta become non-relativistic, i] ~ 350 .See 
table ^ for notation. 

completely different for the case in which the quanta of 
N have a sufficiently high mass. Indeed, when the ampli- 
fication of X from the r.h. neutrino quanta substantially 
decreases when the latter become non-relativistic. As a 
consequence of the two effects mentioned in this para- 
graph, the asymptotic value of Nc.x (at least in the time 
interval we have considered) decreases at small h . Fig- 
ure H^also shows a decrease of the asymptotic Nc^x for 
high h . As discussed in , for sufficiently high coupling 
h the two fields have comparable occupation numbers, 
Nc^N — Nc^x ■ Clearly, the higher the coupling is, the 
sooner this (approximate) equipartition is reached. For 
a very high h , the potential energy associated to the last 
term of (Eq) then disfavors the production of the quanta 
of the two fields, in the same way as a high quartic cou- 
pling cx TV* added to the potential (^ would have pre- 
vented the amplification of the r.h. neutrinos in the two 
fields case. In figure |^, we show the time evolution of the 
effective mass of the quanta of X . Comparing it with 
their distribution in momentum space, one realizes that 
most of the quanta are always in a relativistic regime. 

V. PERTURBATIVE PRODUCTION OF 
GRAVITINOS 

The light quanta X generated at rescattering can in 
turn be responsible for the production of unwanted relics 
such as gravitinos. If unstable, gravitinos with a mass 
of the order of the electroweak scale (which is what we 
expect in models of gravitationally mediated supersym- 
metry breaking) disrupt the successful predictions of pri- 
mordial nucleosynthesis, unless their abundance is below 
the very stringent bound |20[] 

>^3/2 < 10-" (TeV/m3/2) . (29) 
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In inflationary theories, several sources of gravitino pro- 
duction have been considered. The most standard of 
them is the perturbative production from the thermal 
bath formed at reheating. In this case, the above 
limit (|2^) translates into the upper bound Tj-h ^ few x 
10^" GeV on the reheating temperature j2^. Other 
sources of gravitino production can be studied. Since 
gravitinos arise in supersymmetric models, a very natu- 
ral "candidate" channel for their production is the decay 
of the inflaton into its supersymmetric partner (the in- 
flatino) plus a gravitino. This process is however either 
kinematically forbidden ||l^ or strongly suppressed |6|] 
by the fact that the difference between the inflaton and 
the inflatino mass is governed by the supersymmetry 
breaking scale, which is of the same order of the grav- 
itino mass. The resulting gravitino production is suffi- 
ciently small. Recently, the generation of gravitinos at 
preheating has been extensively discussed, both concern- 
ing the relatively easier case of the transverse compo- 
nent , and the more delicate issue of the longitudinal 
component H, 0, ||. These studies are focused on 
the nonperturbative amplification of the gravitino field 
due to the coherent oscillations of the inflaton, and this 
mechanism of production is found to be sufficiently lim- 
ited provided that the inflationary sector of the the- 
ory is weakly (e.g. gravitationally) coupled to the one 
responsible for the present supersymmetry breakdown. 
In this section we discuss a different possible source of 
gravitinos, namely a perturbative production from the 
nonthermal distributions of light MSSM quanta gener- 
ated at rescattering. A comparison with the standard 
thermal production may be used as an initial motiva- 
tion. Concerning the latter, the requirement of a low 
reheating temperature can be seen as the demand that 
the inflaton decays sufficiently late, so that particle in 
the thermal bath have sufficiently low number densities 
and energies when they form. If ~ 10"'^^ GeV at the 
end of inflation, and if the scale factor a is normalized 
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to one at this time, the generation of the thermal bath 
cannot occur before a ~ 10^ . Gravitino overproduction 
is avoided by the fact that in the earHer times most of 
the energy density of the Universe is still stored in the 
coherent inflaton oscillations. As we have seen in the 
previous sections, this last assumption is no longer valid 
if preheating and rescattering effects are important. In- 
deed, in the model considered above the energy density of 
the scalar distributions becomes dominant already when 
the scale factor is of order 100 (the precise value being a 
function of the parameters of the model) . Although 
this comparison is rather suggestive, it is fair to say that 
the computation of gravitino production in the context 
of rescattering is certainly more difficult than the usual 
thermal production. While in the latter case a perturba- 
tive approach can be adopted, and the final result can be 
readily estimated by computing the rate of 2 — s- 2 pro- 
cesses with one gravitino in the final state, rescattering 
is a highly nonlinear phenomenon. In the bosonic case, 
we already remarked that naive perturbative estimates 
poorly reproduce the initial amplification of the fields X . 
Only towards the end of the rescattering stage the num- 
ber densities of the amplified fields become sufficiently 
small so that 2^2 processes become dominant, as the 
(approximate) conservation of the total comoving occu- 
pation number at relatively late times signals |36|| . Un- 
luckily, Pauli blocking forbids fermionic fields to behave 
classically (in the sense discussed in section 3), and lattice 
simulations cannot be used. However, we may take the 
numerical results for bosons as a guideline. Also for the 
production of fermions more complicated processes than 
just 2^2 interactions could be relevant during most of 
the rescattering stage, while they should be subdominant 
at sufficiently late time. The latter is presumably set by 
the same timescale at which rescattering is seen to end 
in the numerical simulations described in the previous 
sections. With this in mind, we proceed to an estimate 
of the amount of gravitinos produced by 2 — s- 2 processes 
with the fields amplified at preheating/rescattering in the 
incoming state. We stress once more that this estimate 
can be reliable only from the beginning of the thermal- 
ization stage on, so it should provide a lower bound on 
the total number of produced gravitinos. It is possible 
that a higher amount of gravitinos is produced at earlier 
times, when nonlinear effects cannot be neglected. This 
computation has been carried out in appendix, where 
also some details (i.e. concerning the quantization of the 
system and the choice to focus on the transverse grav- 
itino component) are reported. As for the thermal pro- 
duction, the dominant 2 — > 2 processes have only one 
gravitino in the final state, and hence only one gravita- 
tionally suppressed vertex. In the thermal case, the dom- 
inant processes have a gauge interaction as the second 
vertex, consider e.g. the process X X ^'03/2 with one 
higgsino in the propagator. For the nonthermal distri- 
butions of scalars that we are considering, however, such 
processes are kinematically forbidden. This is a crucial 
point, which poses a significant limit on the estimated 



production of gravitinos. We can easily understand it 
using the specific process just mentioned as an example: 
either the Higgses are not amplified at rescattering (so 
the above process is irrelevant) or the non- vanishing (i?^) 
provides an effective mass to the zino produced in the in- 
teraction. When the light scalar distributions saturate 
(that is, when the gravitino production can be effective), 

we find numerically a ^/{X^ ^ (lO~^ - 10"^) ■ The 
typical comoving momenta ap characterizing the scalar 
distributions are instead only about one order of magni- 
tude higher than the "inflaton mass" at the end of in- 
flation, ap ^ 15\/A0o • As a consequence, the gauginos 
acquire a mass ^ (10^ — 10'^) p, which shows that 
these processes are kinematically forbidden. We no- 
tice that, at least for this specific kind of interactions, the 
system is still effectively behaving as a condensate: the 
number densities of the scalar distributions are set by the 
quantity (X^) , which is much higher than the typical 
momenta of the distributions themselves. This generates 
a high effective mass for all the particles "strongly" cou- 
pled to these scalar fields. Compare this situation with 
a medium in thermodynamical equilibrium: in this case 
both the typical momenta and the effective masses are set 
by the only energy scale present, namely the temperature 
of the system. As should be clear by the above discussion 
(see also [|6j), the thermalization of the distributions pro- 
duced at rescattering necessarily proceeds through par- 
ticle fusion. Only after a sufficiently prolonged stage of 
thermalization, will the system be sufficiently close to 
thermodynamical equilibrium so to render processes as 
the one discussed above kinematically allowed. As we 
discuss in the appendix, kinematically allowed 2^2 
interactions can be obtained by taking a trilinear inter- 
action coming from the superpotential term (p6|), also 
responsible for the Dirac mass term for the neutrinos. 
This can lead to processes of the kind NrX xip^/2^ 
or X X ^ Nfi'ip^/2- The physical number density of 
(transverse) gravitinos produced by these processes can 
be estimated as 



^^3/2 iv) 



Ml 



a iv*) 



T] > V* 



where rj^, is the time (to be determined numerically) 
at which the expression in the first parenthesis reaches 
its maximum, while the second parenthesis is a dilu- 
tion factor due to the expansion of the Universe at later 
times (see the appendix for details). As we have re- 
marked, this result is subject to the limit (p9|), where 
= J ^J^d s is the entropy density of the Uni- 

verse, computed once the dominating thermal bath is 
formed. For practical use, we find that a more "conve- 
nient" bound can be obtained if (29) is combined with 
the result for the baryon asymmetry, eqn. (p3|). For this 
purpose, consider the ratio 



c 



^3/2 _ 23 N3/2 iv) _ 23 jV3/2 ij]*) 
Yb ~ 8ei iV^v (77) ~ 8ei TVjv (?]*) 



(31) 
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FIG. 6: Contour plot for the quantity (/h , where ( is defined 



m eqn. 



The contour fines range from (^/h = 3-10 



to 



C//i = 1.1 • 10"^° (darker region). 



The quantity ^ has two main advantages, (i) it is inde- 
pendent of the entropy of the Universe and (ii) it can be 
comptited already at = 77* , since after this time the two 
physical number densities ^^^d simply rescale as 
. While Y3/2 must satisfy the upper bound ([29|), the 
limit Yb ft 10^^^ poses a lower bound on the number 
density of the sneutrinos, if leptogenesis is assumed to be 
responsible for the baryon asymmetry of the Universe. 
Adopting the parameterization we then see that the 
ratio C has to satisfy 



3 X 10*^ 



'JCP 



10^° GeV \ 
M J 



0.05 eV\ A^3/2(?7* 



m3 J NNiv*) 



< 10' 



(32) 



It is important to stress that, unlike the limit (29) 



this bound cannot be ameliorated by an eventual en- 
tropy release which may occur between the decay of the 
r.h. sneutrinos and nucleosynthesis, since both the grav- 
itino and the baryon number densities would be diluted 
in the same amount. For this reason, we find in the 
present context the bound ( p2| ) more significant than 
the limit (29) involving the gravitino abundance alone. 
We wanted to verify whether the condition (^) is re- 
spected for the choice M = 10^^ GeV considered in the 
previous sections, and for several values of the couplings 
g and h (defined in the potential (p8|)) in the range 
g e [30, 5000] , h e [2000, 200000] (we remind that 
g = g^/X, and analogously for h). To do so, we have 
defined 



^ 1 /a^ev^ . 

OCP V ni3 



(33) 



and in figure ^ we have plotted the quantity (/h. In 
this way, we factor out the explicit dependence of C on 
h coming from the cross section of the dominant 2^2 
processes, see e.g. eqn. (^o|). The qualitative behavior of 
with h (vertical axis) is easily understood in terms 
of the arguments used to explain the results of figure ||. 
For relatively low h , the amplification of the X field is 
weak, and so few gravitinos quanta are produced. As 
h increases, the amplification becomes stronger, both re- 
garding the final value at which Nx saturates and the ra- 
pidity at which the saturation occurs. As a consequence, 
the number density of produced gravitinos also increases. 
As h further increases, the rapidity at which X saturates 
keeps increasing (and so, the time rj^ at which most of 
the gravitino quanta are produced decreases); however, 
the final value of Nx starts to become smaller, leading 
to the decrease oi (/h that we observe in the figure. For 
fixed values of h, the final result is also first increasing 
and then decreasing with g . This behavior is presum- 
ably due to the dependence on g of the total number 
of quanta of both the N and X fields produced at pre- 
heating/rescattering (notice that it vanishes both at very 
small and very high g). However, the interpretation is 
in this case less clear. Concerning the value of C itself, 
in the range of coupling considered it reaches the maxi- 
mal value at g ~ 100 , h ~ 200000 , where it evaluates to 
^ ~ 10~^ . From what we have just said, we expect that 
C starts decreasing at higher h , although the numerical 
simulations we have performed show that the decrease 
starts only at the highest value of h that we have con- 
sidered. From the definition (^3|), we see that the bound 
C < 10"'^ is respected, provided the CP violation en- 
coded in the parameter Sqp (see eqs. (^) and (||)) is not 
too small. However, we remind that our estimates take 
into account only the gravitino quanta produced from the 
end of the rcscattering stage on, while a higher produc- 
tion at earlier times cannot be excluded. 
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APPENDIX 

In this appendix we derive eqn. (|30| ) of the main text. 
We first estimate the number density of gravitinos pro- 
duced by the nonthermal distributions formed at rescat- 
tering. We remark that at this stage supersymmetry 
breaking is controlled by the energy density of these dis- 
tributions. The longitudinal gravitino component (i.e. 
the goldstino, which is "eaten" in the unitary gauge) is 
thus provided by a linear combination of the fermionic su- 
perpartners of these scalars, and does not coincide with 
the longitudinal gravitino at late times (at least for the 
standard case of a present gravitationally mediated su- 
persymmetry breakdown). For this reason, our discus- 
sion will be limited to the transverse gravitino component 
"03/2 • Its mass is given by 

-3/2 = e-(«A4^^, (A-1) 

where, we remind, W and /C are, respectively, the super- 
potential and the Kahler potential of the model, while 
denotes the set of scalar fields amplified during pre- 
heating and rescattering. To quantize the transverse 
gravitino, we define a homogeneous mass rriy^ by 
placing in eqn. ( |A-l| ) the (x— dependent) values of the 
fields with their (homogeneous) variances, (/)(t, x) 
^ (t) . The difference bm-j^ji = m-j,i2 — be 
accounted for in the interaction lagrangian. The main 
production of gravitinos is expected to occur close to 
the point at which the number density of light scalars 
X reaches its maximum, in the same way as most of 
the thermal production occurs as soon as the thermal 
bath is generated. From the results of section 4, we ob- 
serve that the maximum of Nc^x is achieved at the end 
of the rescattering phase. At this moment the thermal- 
ization stage begins, during which the variances of the 
fields show an adiabatic evolution. This allows us a con- 
sistent quantization of the gravitino component, since the 
mass 7713^2 {^) is also varying only adiabatically in this pe- 
riod. It is clear that the main concern with the procedure 
just described is that, contrary to the thermal case, the 
difference 5m^/2 is now of the same order of m^^r^ it- 
self, at least during the initial part of the thermalization 
stage. This leads to the problem discussed in section 5, 
namely to the fact that the perturbative computation of 
the gravitino production is presumably not as straight- 
forward as in the thermal case, and more complicated 
processes than just 2 — > 2 interactions can be expected 
to be relevant. However, as we have already remarked, 



the latter should provide at least an order of magnitude 
estimate for the gravitino produced from the end of the 
rescattering stage on, and should reasonably lead to a 
lower bound to the total production. For this reason, we 
now proceed to a rough estimate of their cross sections. 
The dominant processes with two gravitinos as outgoing 
particles have two gravitationally suppressed vertices (i.e. 
XX —> "03/2 "03/2 with a flat direction fermion x in the 
propagator; processes coming from the interaction term 
1^^3/2 03/2 03/2 are subdominant). Their cross section is 
of the order a ^ /Mp , where here and in the following 
p denotes the typical momentum exchanged in the scat- 
tering. As in the thermal case, the distributions of the 
light quanta are indeed characterized by a typical mo- 
mentum; while for a thermal distribution p ~ T , we now 
have p ~ K VA 0o (t) 7 where in the cases shown below 
K is a coefficient of order 10 dependent on the specific 
choice of the parameters. In our estimates we will take 
K ~ 15 . @ Thus, cr ~ 10"" Mp^ for this class of 
processes. A more efficient production is expected from 
scatterings with only one gravitationally suppressed ver- 
tex, and hence with only one gravitino in the final state. 
For example, the standard thermal production is mainly 
due to channels having a gauge interaction as the second 
vertex, e.g. H H ^ 03/2^ with an exchanged higgsino. 
In the present context, however, processes with outgoing 
gauginos (that we generically denote with g) are expected 
to be kinematically forbidden, since these particles ac- 
quire a high effective mass from their interaction with 
the nonthermal scalar distributions. Indeed, if a scalar 
field X has a large vev, and an interaction of the form 
y/aX^g is present (0 and g are two component matter 
fermion and gaugino, we use ^/a because we have already 
used g as the inflaton-neutrino coupling) then the gaug- 
ino acquires a large Dirac mass ~ \/oi{X) mixing with 0. 
We have large variances, rather than a large vev; by anal- 
ogy with finite temperature, we expect that (X^) ^ will 
generate an effective mass square "m^ ~ a(X^)" in the 
g and dispersion relations. Q So for kinematic pur- 
poses, we assume that gauginos which couple to the fiat 
direction have masses of order ^Ja{X'^). When the light 
scalar distributions saturate (that is, when the gravitino 
production we are discussing can be effective), we find 
numerically (X^) ^ (lO^^ — 10~^) <f)o/a . As a conse- 
quence, rrig ~ (10^ — 10'^) p, and these scatterings are 
forbidden. One is immediately led to consider processes 
with an additional Xi ipj g vertex and in which the heavy 
gaugino is off-shell. Their cross section can be roughly 
estimated as cr 10-2 {a/Mpf {p/mgY , which IS com- 
parable or smaller than the cross section for the process 
XX —> 03/2 03/2 considered above. Finally, there is the 
possibility that the second vertex comes from the super- 
potential term (p6|), also responsible for the Dirac mass 
term for the neutrinos. This can lead to processes of 
the kind NrX ^ a; 03/2 or XX ^ -^iJ0'3/2 {x denot- 
ing the fermionic partner of X ; all processes have in the 
propagator the fermionic partner of one of the incom- 
ing scalars). The cross sections for these processes are 
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roughly estimated as ||§ cr - /i^/Mj, ~ IQ-^^ /i/M|,. 
Thus, unless of a very small coupling h , the last class of 
scatterings has the highest cross section and dominates 
the production of the transverse gravitinos. In particu- 
lar, processes with one incoming Nn quantum are dom- 
inant if Nc^x starts to saturate at a smaller value than 
Nc.N ■ Viceversa, scatterings of the kind X X Nu -03/2 
will dominate. Numerical results indicate that the for- 
mer situation is more often realized. The cases in which 
the opposite was found are characterized by a relatively 
high coupling h , so that the light degrees of freedom are 
quickly amplified to values Nc,x ^ -^c.at ■ Anyhow, in 
these cases the cross sections of the two type of processes 
are clearly of the same order of magnitude. Hence, for 
brevity of exposition we will only refer to the processes 
with an incoming r.h. sneutrino, although both the two 
possibilities have been considered in our estimates. The 
integrated Boltzmann equation reads 



3/2 



dt 



+ 3HN3f2-{'T\v\)NxNN , 



(A-2) 



where the "friction term" due to the expansion of the 
Universe can be neglected in the estimate of the order of 



magnitude of gravitinos produced. The whole production 
time can be then divided in a series of time intervals of 
duration H~ ^ (ti) each. During each interval, quanta of 
gravitinos are generated with a density of 



a{ri^) a{rii) 



(A-3) 



(notice the presence of the scale factor, since the physi- 
cal and not the comoving occupation number has to be 
used in the integrated Boltzmann equation). The func- 
tion i^Ncx Nc^N a~^) (77) amounts to zero at the end 
of inflation, and it reaches a maximum at a time 77* , 
which can be determined numerically and which roughly 
corresponds to the moment at which the comoving num- 
ber density Nc^x starts saturating (this in turns occurs 
towards the end of the rescattering stage, when 2^2 
processes start dominating). At > 77* it then quickly 
decreases due to the expansion of the Universe. It thus 
turns out that, as for the thermal production, the grav- 
itino quanta are mostly generated at the time 77* , so that 
their "late time" physical number density is approxima- 
tively given by 
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N3/2 iv) 



Ml 



a (77* 



a (77*) 
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n 3 



a (77) 



(A-4) 



This is eqn. ( ^ ) of the main text. 
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We will use A'^ as a shorthand for the superfield, its scalar 
and fermionic component. We explicitly refer to the par- 
ticle type when this could cause a confusion. 
A nonperturbative inflaton decay also occurs for hybrid 
inflation [^, However, we will not discuss this 

possibility here. 

In some works, the term rescattering refers only to the 
scatterings of the produced quanta against the zero mode 
of the inflaton. Here we keep the original meaning given 



in [0. 

[77] Gauge fields which are not coupled to the fields generated 
at rescattering will not acquire this high effective mass. 
However, being uncoupled, they will not be amplified ei- 
ther. 

[78] The situation is even more enhanced for hybrid inflation- 
ary model, in which the energy density of the zero modes 
of the scala rs g ets dissipated within their first oscilla- 
tion [||, H |of . 

[79] We acknowledge very useful discussions with Patrick B. 
Greene and Lev Kofman on this issue. 

[80] The superfield A*' written in eqn. (1) actually denotes an 
anti-(s)neutrino. However, for brevity we will refer to it 
as a (s)neutrino. 

[81] Any subsequent entropy production leads to further di- 
lution of the asymmetry. 

[82] We will use the parametrisation (^) for all the ei, i = 1..3. 
It is possible that Sep < 1 for 62 and es (assuming no 
cancellations in the formulae), although this has not been 
shown. 

[83] This estimate is not significantly changed if the angles 
in Vl are small compared to ^13. CHOOZ experiment: 

013 .1. If 013 < .1, but 013 < [VlU2,[Vl]i3 < -1, 

then the formula for e is similar to WW, with the re- 
placement 6ii Oli] and S (pij {where [Vl]i2 ~ 
cos 01,13 sin 6•il2e'»'l^ [Vl]i3 ^ smOLise"^'^). 

[84] The "almost" is important; if W — 1, there is no CP 
violation, so e = 0. 

[85] We assume in this work that Trh is "large"; for a dis- 
cussion of baryogenesis in Iow-Trj/ models, see e.g. [^ . 

[86] One should also consider the perturbative decay of the 
inflaton quanta. Comparing eqn. (^6|) with the number 
of r.h. neutrinos produced perturbatively in one inflaton 
oscillation (i.e. the typical timescale for preheating), one 
can however see that the perturbative production is sub- 
dominant when kinematically allowed. 

[87] For instance, this is the case if the r.h. neutrinos lifetime 
is long enough for them to come to dominate the energy 
density in the Universe, after reheating has completed. 

[88] The thermalization of this system is a very interesting 
issue, which however we do not discuss in this work - 
see [Q for a more detailed study. Since in this case 
thermalization proceeds via particle fusion, an impor- 
tant role may be played by three or five point vertices, 
which shorten perturbative estimates of the thermalisa- 
tion timescale |60|] (other recent discussions on thermal- 
ization can be found in |^). 

[89] To embed the system in a supergravity context while pre- 
serving a flat potential for the inflaton field, one may 
impose [^ a definite parity for the Kahler potential 
= ("I> — "!>*) . Doing so, the inflaton is identified with 
the real direction of the scalar component of $ , and su- 
pergravity corrections can be neglected. 

[90] We neglect the term quartic in A'' , subdominant with 
respect to the mass term N"^ during most of the pre- 
heating/rescattering period, as well as the mixed term 
oc \f\ g (fy^ , negligible with respect to the mixed term 
present in eqn. (K% for 2> A . One may also be worried 
that, if the right-handed sneutrino is charged under some 
gauge group (as it generally happens in grand unified 
models) with a gauge coupling go not much smaller than 
one, the corresponding D— term oc gQ \N\'^ C V could 
prevent the amplification of A*' at preheating-rescattering. 
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However, at least as long as (A'^) is smaller than the scale 
at which the gauge symmetry is broken, this term gets 
actually compensated by a shift of the (much heavier) 
field responsible for the breaking of the symmetry, and it 
is thus absent from the effective potential for A'' Q . 
[91] Non-adiabatic production of sneutrinos can occur for 
bare masses as large as p 0o/4 ~ g ■ 10^* GeV, but the 
efficiency of the process will be much lower, because red- 
shift effects will terminate the resonance before rescat- 
tering sets in. 

[92] The precise value of the typical momentum fe« of the 
distributions, as well as the time needed for Nc to satu- 
rate, are a nontrivial function of g , since different values 
of this parameter lead to different positions (in momen- 
tum space) and strengths of the resonance bands [pif . 
However, the rescattering stage destroys these resonance 
bands, making the dependence of fe* on g milder. 

[93] Notice that for a massless inflaton, the inflaton energy 
redshifts as radiation, and non-relativistic neutrinos will 
easily dominate the energy in the Universe. If the inflaton 
in instead massive, then the energy in sneutrinos would in 
any case be a fraction of order unity of the background 
energy, and would start increasing after inflaton decay. 
As a consequence, the resulting baryon asymmetry would 
still be givenfat most up to factors of order one) by the 
expression (Ea). 

[94] The situation is even more enhanced for hybrid inflation- 
ary model, in which the energy density of the zero modes 



of the scala rs g ets dissipated within their first oscilla- 
tion [||, H |of. 

[95] Processes with an additional gauge interaction and the 
gaugino in the propagator are allowed but strongly sup- 
pressed. See the appendix for details. 

[96] The time rj, roughly corresponds at the moment at which 
the distribution of light quanta X starts to saturate. This 
typically occurs towards the end of the rescattering stage, 
which guarantees that considering only 2 — > 2 processes 
should provide at least an order of magnitude estimate 
of the gravitino quanta produced at this stage. 

[97] The existence of a typical momentum allows the use 
of the integrated Boltzmann equation to estimate the 
amount of gravitinos produced. Moreover, since this mo- 
mentum is much higher than the gravitino mass, the 
value of the latter does not affect the cross sections for 
the processes with outgoing transverse gravitinos. This is 
welcome, since the above (somewhat arbitrary) redefini- 
tion m3/2 ~ rn%/2 + ^'^3/2 will not affect our estimates. 

[98] ft is implicit in this discussion that all the fermionic fields 
are quantized in the same way as done for the gravitinos. 

[99] In this estimate it is assumed that the exchanged momen- 
tum p is higher or at most comparable with the mass of 
the r.h. neutrinos. This is certainly true when most of 
the gravitinos are produced, i.e. when the distributions 
of light quanta X are about to saturate. 



